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b 1 Abstract 
< 

. The Schrodinger functional formalism is given as a field theory in a finite 

! volume with a Dirichlet boundary condition in temporal direction. When 

£> \ one tries to construct this formalism with the Ginsparg- Wilson fermion 

^ ; including the overlap Dirac operator and the domain-wall fermion one 

O ' easily runs into difficulties. The reason is that if the Dirichlet boundary 

condition is simply imposed on the Wilson Dirac operator Dw inside 
! of the overlap Dirac operator an exponentially small eigenvalue appears 

in Dw, which affects the locality properties of the operator. In this 
paper we propose a new procedure to impose the Schrodinger functional 
Dirichlet boundary condition on the overlap Dirac operator using an 
orbifolding projection. 
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I. INTRODUCTION 



The Schrddinger functional (SF) is defined as a transition amplitude between two bound- 
ary states with finite time separation [1,2] 

Z = (C; x = T\C; x = 0) = J VQe'^ (1.1) 

and is written in a path integral representation of the field theory with some specific bound- 
ary condition. 1 One of applications of the SF is to define a renormalization scheme beyond 
perturbation theory, where the renormalization scale is given by a finite volume TxL 3 ~ L 4 
of the system. The formulation is already accomplished for the non-linear a- model [5], the 
non-Abelian gauge theory [6] and the QCD with the Wilson fermion [7,8] including 0(a) 
improvement procedure [9,10]. (See Ref. [11] for review.) 

In this formalism several renormalization quantities like running gauge coupling [12-17], 
Z-factors and 0(a) improvement factors [18-22] are extracted conveniently by using a Dirich- 
let boundary conditions for spatial component of the gauge field 

A k (x)\ Xo=0 = C k (x), A k (x)\ X0=T = C' k (x) (1.2) 

and for the quark fields 

P+ip(x)\ X0=0 = p(x), P_ip(x)\ X0=T = p'(x), (1.3) 
ip(x)P-\ X0=0 = p(x), i/j(x)P+\ X0=t = p'(x), (1.4) 

P± = ^ d.5) 

One of privilege of this Dirichlet boundary condition is that the system acquire a mass gap 
and there is no infra-red divergence problem. 

We notice that the boundary condition is not free to set since it generally breaks symme- 
try of the theory and may affect renormalizability. However the field theory with Dirichlet 
boundary condition is shown to be renormalizable for the pure gauge theory [6]. And it is 
also the case for the Wilson fermion [8] by including a shift in the boundary fields. 

Although it is essential to adopt Dirichlet boundary condition for a mass gap and renor- 
malizability, it has a potential problem of zero mode in fermion system. For instance starting 
from a free Lagrangian 

£ = ^M + m)ij (1.6) 
with positive mass m > and the Dirichlet boundary condition 

P_^U=o = 0, P + VL=t = (1.7) 
the zero eigenvalue equation (7o<9q + m)ip = in temporal direction allows a solution 



1 The same kind of transition amplitude was introduced in [3,4] in order to implement the temporal 
gauge. 
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$ = P +e - mx ° + p_ e -™(T-xo) ( L g) 

in T — > oo limit and a similar solution remains even for finite T with an exponentially small 
eigenvalue oc e~ mT . In the SF formalism this solution is forbidden by adopting an "opposite" 
Dirichlet boundary condition (1.3) and the system has a finite gap even for m = [7]. 

In the SF formalism of the Wilson fermion [7] we cut the Wilson Dirac operator at the 
boundary and the Dirichlet boundary condition is automatically chosen among 

P±i>\x =o = 0, P^U=t = (1.9) 

depending on signature of the Wilson term. For example if we adopt a typical signature of 
the Wilson term 

d w = i,\ (v; + v,) - ^v;v, + m (i.io) 

the allowed Dirichlet boundary condition is the same as (1.3). In this case the zero mode 
solution can be forbidden by choosing a proper signature for the mass term; the mass should 
be kept positive M > to eliminate the zero mode [7]. 

However this zero mode problem may become fatal in the Ginsparg- Wilson fermion 
including the overlap Dirac operator [23,24] and the domain- wall fermion [25-28]. The 
overlap Dirac operator is defined by using the Wilson Dirac operator (1.10) as 



! > =M D„ -p=_ |. a=— . (1.11) 

Here we notice that the Wilson fermion mass M should be kept negative in a range —2 < 
M < to impose heavy masses on the doubler modes and a single massless mode to survive. 
As explained in the above when the Dirichlet boundary condition is imposed directly to 
the kernel Dw an exponentially small eigenvalue is allowed for this choice of the Wilson 
parameter and the mass. Nearly zero eigenvalue in Dw may break the locality of the overlap 
Dirac operator [29]. The situation is also quite similar in the domain- wall fermion. The same 
zero mode solution appears in the transfer matrix in fifth direction, which suppresses the 
dumping solution in fifth dimension and allow a chiral symmetry breaking term to appear 
in the Ward-Takahashi identity [27]. 

Since a naive formulation of the SF formalism by setting the Dirichlet boundary condition 
for the kernel Dw does not work, we need different procedure to impose boundary condition 
on the overlap Dirac fermion. In this paper we propose an orbifolding projection for this 
purpose. In section II we introduce an orbifolding in a continuum theory and show that the 
Dirac operator of the orbifolded theory satisfy the same SF boundary condition. Although 
this procedure may be applicable to the general Ginsparg- Wilson fermion we concentrate on 
the overlap Dirac fermion in section III. Section IV is devoted for conclusion. 



II. ORBIFOLDING FOR CONTINUUM THEORY 

We notice a fact that the chiral symmetry is broken explicitly by the Dirichlet boundary 
condition (1.3) in the SF formalism. This should be also true in the overlap Dirac operator; 
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the Ginsparg- Wilson relation should be broken in some sense, which was not accomplished 
in a naive formulation. We would adopt this property as a criterion of the SF formalism. 

Then we remind a fact that an orbifolded field theory is equivalent to a field theory with 
some specific boundary condition. Since it is possible to break chiral symmetry by an orb- 
ifolding projection in general, it may be able to represent the SF formalism as an orbifolded 
theory. In this section we search for an orbifolding projection which is not consistent with 
chiral symmetry and provide the same SF Dirichlet boundary condition (1.3) and (1.4) at 
fixed points. 

We consider a massless free fermion on S 1 x R 3 

C = ^(x)^d^(x), (II.l) 
where the anti-periodic boundary condition is set in temporal direction of length 2T 

ij;(x, x + 2T) = -tp(x, x ), ij;(x, x + 2T) = -tp(x, x ). (11.2) 

The orbifolding S 1 /Z 2 in temporal direction is accomplished by identifying the negative 
time with the positive one x <-> — x . Identification of the fermion field is given by using a 
symmetry transformation including the time reflection 

il>(x) -> T,il>(x), if;(x) -> i/>(x)E, E = i7 5 7o-R, (H.3) 

where R is a time reflection operator 

Rip(x, x ) = if>(x, -x Q ). (II.4) 

R has two fixed points Xq = 0,T, where Xq = is a symmetric and xq = T is an anti- 
symmetric fixed point because of the anti-periodicity 

Rip(x, 0) = ip(x, 0), Ri/j(x, T) = -rj}{x, T). (11.5) 

It is free to add any internal symmetry transformation for the identification and we use the 
chiral symmetry of the massless fermion 

ip(x) -> -ij 5 tp(x), ip(x) -> -ip(x)i^ 5 . (11.6) 

Combining (II. 3) and (II. 6) we have the orbifolding symmetry transformation 

ijj(x) -> -Ity(ar), ^(x) -> ^(x)T, T = 70J R. (II.7) 

The orbifolding of the fermion field is given by selecting the following symmetric sub- 
space 

n + ^(x) = o, ^(x)n_ = o, n ± = i^. (ii.s) 

We notice that this orbifolding projection provides the proper homogeneous SF Dirichlet 
boundary condition at fixed points Xq = 0,T 
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P + ^(x)| Xo=0 = 0, P_^(x)| Xo=T = 0, (II.9) 

iP(x)P_\ X0=0 = 0, 4>(x)P + \ X0=T = 0. (11.10) 

The orbifolded action is given by the same projection 

s = ±J d 4 x^(x)D SF ij(x), d sf = n + ^n_, (ii.il) 



where factor 1/2 is included since the temporal direction is doubled compared to the original 
SF formalism. We notice that the chiral symmetry is broken explicitly for Dsf by the 
projection. 

Now we have two comments. Since the Schrodinger functional of the pure gauge theory 
is already well established [6] we treat the gauge field as an external field and adopt a 
configuration which is time reflection invariant 

A (x, -x ) = -Ao(x,x ), Ai(x, -x ) = Ai(x,x ) (H.12) 

and satisfy the SF boundary condition (1.2) simultaneously. We set periodic boundary 
condition for the gauge field 

A fl (x,x + 2T) =A lx (x,x ). (11.13) 

Second comment is on the mass term. Although the mass term is not consistent with the 
chiral symmetry, we can find a symmetric mass term under the orbifolding transformation. 
A requirement is that the mass matrix M should anti-commute with orbifolding operator 
{M, T} = 0. One of the candidate is a time dependent mass M = mr](xo) with anti- 
symmetric and periodic step function 

rj(-xo) = -ri(xo), i](x + 2T) = r](x ), 

r](x ) = 1 for < x < T. (11.14) 
Now the Dirac operator becomes 

Dim) = 7 M (<9 M - iA^x)) + mrj(x ), 
which has the orbifolding symmetry 

{D(m),T} = 

The orbifolded action is defined in the same manner as the free theory 

S = i y d 4 x$(x)D SF (m)ij(x), D SF {m) = U + D{m)U_. (11.17) 

We can show that this system with gauge interaction is equivalent to the original QCD 
with SF boundary condition [7]. We consider the bulk region of the temporal direction: 
— T < x < and < x < T. In this region the orbifolding condition (II. 8) becomes 

ip(x ) = --fotp(-xo), ip(x ) = ^(-xoho (11.18) 
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to identify the fields in positive and negative time. Combining with the anti-commutative 
property of the Dirac operator (11.16) we can show that contribution to the action from the 
bulk region is written as 

Sbuik = \ I d 4 x^(x)U + D(m)U^(x) 

2 J bulk 

= f T dx d 3 xiJ(x) [ 7/i (<9 m - iA^x)) + m\ ij)(x), (11.19) 

which is nothing but the QCD action introduced in Ref. [7] except for the surface term, 
which will be discussed later. Since the bulk action and the boundary condition for the 
orbifold construction is exactly the same as those of QCD with SF boundary condition, the 
spectrum and the propagator should be uniquely determined to be equivalent to those of 
Ref. [7] and Ref. [10]. One can easily check this fact at tree level. 

III. ORBIFOLDING FOR OVERLAP DIRAC FERMION 

Application of the orbifolding procedure is straightforward to the Ginsparg- Wilson 
fermions including the overlap Dirac operator [23,24], the domain- wall fermion [25-28] and 
the perfect action [30-32] which possess both the time reflection symmetry 

[£,D] = (III.l) 

and the lattice chiral symmetry [33] stemming from the Ginsparg- Wilson relation [34] 

7s D + D l5 = aD l5 D. (III.2) 

In this subsection we concentrate on the overlap Dirac operator (1. 11), for which the time 
reflection symmetry (III.l) comes from that of the Wilson Dirac operator [E, Dw] = 0. 

A. Orbifolding construction of Dirichlet boundary 

As in the continuum case we consider a massless fermion on a lattice 2N T x JV| with 
anti-periodic boundary condition in temporal direction (II. 2). We use an orbifolding S 1 /Z 2 
in temporal direction. Identification of the fermion field is given by using the time reflection 
(II. 3) and the chiral symmetry of the overlap Dirac fermion [33] 

ip(x) -> -i%ip(x), ip(x) -> -^(x)i7 5 , 75 = 75 (1 - 3D) , (III.3) 

where the gauge field is treated as an external field and we adopt a time reflection symmetric 
configuration 

U k (x,x ) = U k (x, -x ), U (x,x ) = Ul(x, -x Q - 1), (III.4) 
satisfying the SF Dirichlet boundary condition simultaneously 

U k (x, 0) = W k (x), U k (x, N T ) = W' k {x). (III.5) 
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Combining (II. 3) and (III. 3) we have the orbifolding symmetry transformation 

ip(x) -> -TiP(x), ${x) -> lj>(x)r, f = T(l - aD), (III.6) 

where T is the same as the continuum one (II. 7). We notice that starting from the time 
reflection symmetry of the Dirac operator (III. 1) and the Ginsparg- Wilson relation (III. 2) 
we have another GW relation for T 

YD + DY = aDYD (III.7) 

and T Hermiticity 

YDY = DK (III.8) 

The operator Y has a property Y 2 = 1 like Y and can be used to define a projection operator 
in the following. 

The orbifolding identification of the fermion field is given in the same way with slightly 
different projection operator 

n+^(x)=o, ^(z)n_ = o, n± = ^^, (111.9) 

which turn out to be the SF Dirichlet boundary condition (II. 9) and (11.10) at fixed points 
in the continuum limit. Using the time reflection symmetry (III. 1 ) we can easily show that 
the projection operators Y and Y do not have an "index" 2 

t r r = trf = (III. 10) 

and furthermore we can find a local unitary transformation 

u = I±^(l-aZ>) + I^, u' = l5 u l5 , (111.11) 
which connects Y and Y as 

f = «tr«, f = u'Yu'l (111.12) 

The projection operator II ± spans essentially the same Hilbert sub-space as IT ± . We notice 
that this unitary operator connects % and 75 in a similar way 

75 = u^ 5 u, 75 = u'-f 5 u n . (III. 13) 

The physical quark operator is defined to transform in a same manner as the continuum 
under chiral rotation, 

Sq(x) = 75<?(x), 5q(x) =q(x)'j 5 . (III. 14) 



2 The same property is satisfied for the chiral index tr% = tr^E 2 ) = — tr75 = 0. 
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Since we have a unitary operator u and u' we have several ways to define a physical quark 
field from GW fermion fields ip and ijj. 3 For example 

q(x) = (l- |l>) 4>(x), q(x) = #r), (111.15) 

= uip(x), q(x) = ip(x), (III. 16) 

q(x) = u'^ip(x), q(x) = ip(x). (III. 17) 

These three definitions are not independent but connected with u + u'^ = (2 — aD). The 
orbifolding of the physical quark field becomes the same as that of the continuum theory 

Il + q(x) = 0, g(x)n_ = 0. (111.18) 

The massless orbifolded action is given by 

S = ^ 4 E^sf^, -Dsf = H+DT1-. (111.19) 

We have four comments here, (i) It should be emphasized that the SF Dirac operator D SF 
is local since it is constructed by multiplying local objects only, (ii) The massless SF Dirac 
operator D$f does not satisfy the chiral Ginsparg- Wilson relation (III. 2) and the chiral 
symmetry is broken by projection as was expected, (iii) Although two different projection 
operators T and T are used from the left and right of D$p this does not bring the problem 
we encountered in the chiral gauge theory since these two operators are connected by the 
unitary transformation u or u' . (iv) Since we have some ambiguity in defining the mass 
term we concentrate on the massless case and postpone the massive theory to the later 
sub-section. 

It is clear that (III. 19) is a well defined lattice regularization of the SF QCD with gauge 
interaction since the Dirac operator is local and the lattice action (III. 19) converges to the 
continuum one (11.17) in a — > limit including the boundary condition. We shall see explicit 
scaling behavior for free theory in the following two sub-sections. 

B. Eigenvalue of free Dirac operator 

We consider an eigenvalue problem of the massless SF Dirac operator 

£> S f = n+DlL (111.20) 



3 One may suppose that a local Dirac operator can be constructed with exact chiral symmetry 
by adopting the physical quark operator of (III. 16) and (III. 17), where the effective action of q, q 
is kept local. However the unitary transformation u contains the time reflection operator R and 
furthermore the gauge configuration is fixed to be time reflection invariant. There is no translation 
invariance in the physical quark effective theory and this is consistent with the Nielsen-Ninomiya 
no-go theorem. 
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at tree level. This Dirac operator connects two different Hilbert sub-space as in the contin- 
uum [7] 

D SF : H + - H-, D\ ¥ : - H+, (111.21) 

ft + = {?/>|n + ?/> = o}, H- = {^|n_^ = 0}. (111.22) 

We need to introduce the "doubled" Hermitian Dirac operator 

*=U DsT )' (IIL23) 

which connects the same Hilbert space H- © H + — > 7i_ © 7i + . As in the continuum the 
eigenvalue problem is solved on a two component vector 

V = (t-), ^fj.eH-, $+eH+ (111.24) 



and the eigenvalue equation is written in a following form 

£>sf^ + = A^_, 4f^- = (111.25) 

with a real eigenvalue A. 

A candidate of the eigen-function in each Hilbert sub-space H + and is given by 

Mx) = uH + {x), jP + (x) = f + (x )ve ifB , = f_{x Q )wj&, (111.26) 

where /± are given to satisfy r/± = =f/± as 

f+(x ) = a (P+ sinpo^o + P- cosp ^o) , (III. 27) 

f-(x ) = a (P + cosp x + P- sinp x ) (III. 28) 

with 

2n — 1 

Po = ~2NjT n ' n = i- N T + ^---iN T } (HI.29) 

for anti-periodicity in 2A^. t> and w are four component vectors in spinor space. In the 
following we consider the operator 

zW = n + {±(p-D*)- t}- (d + £>t) I n_ (111.30) 

and its operation on 

We start from the fact that ip + is an eigen-function of the operator D\yD w and D w + D\y 
with eigenvalues A 2 ^ and 2W where 

A^ = £sin 2 ap M + W 2 , W = -\M\ + £ (1 - cos ap M ) . (111.31) 
On the other hand D w — operates as 
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^ (p w - D\^j ip + {x) = (smap f-(x ) + hi sin apif+(x )) ve ipx . (III. 32) 

Combining these two and taking similar manipulation for uDl F on ip_ the eigenvalue equa- 
tion is re- written as 

D SF $ + (x) = /_(x ) (A + inA, + t l5 B) ve^ = X^-(x), (111.33) 

4 F V>- (x) = u f /+ (xo) (A - i 7i A - ^75^) «;e ipB = (x) (111.34) 

with 

Here we used E = i^T and a property r/± = Tf±- At last we solve the eigenvalue equation 
in the spinor space and get the eigenvalue 

A 2 = A\ + B 2 . (111.36) 

This eigenvalue agrees with that of the massless continuum theory [7] in a — > limit com- 
bined with the discretization condition (III. 29). 

C. Free propagator 

The fermion propagator is formally given by the inverse of the Dirac operator 

G SF (:r, y) = 2 (D£) = 2 (fiin + ) . (HI.37) 

x >y \ L) J x ,y 

where inverse is defined in the Hilbert sub-space H + or H- as 

D SF Ds£ = n+, D^Dsp = n_. (IH.38) 
At tree level this propagator can be written in a simple form as 

G SF = 2£>t n+ ^L_n + = £)t (P + G L + P_G R ) , (IH.39) 

G R/L (x,y) = -L E^ ( ^G R/L (x ,y ;p1, (IH.40) 



iV3 

p 

1 I 

G R (x , y ; p) = — — J; inm^ (e* ^^ - e*»<*°-*»>) , (111.41) 

i 1 

G L (x , y ;p) = ^ E 7^—^ (e*>(*°-*»> + e *><*°+«»>) , (111.42) 

^ ?i=— A^+l A 1 

where p satisfy the quantization condition (III. 29). 

If we adopt the definition (III. 15) for the physical quark field its propagator is given by 
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(q(x)q(y)) 



?£> ) G SF 



(x,y). (111.43) 



2 

The proper Dirichlet boundary condition [10] is trivially satisfied for this quark propagator 

P + {q{x)q(y)) U =o = 0, P_ (q(x)q(y)) \ XQ=T = 0, (111.44) 
{Q(x)Q(y)) L=oP- = 0, (q(x)q(y)) L=tP + = (111.45) 
because of the projection IT-|-. At tree level the propagator takes the form 

xe i PO x | (e-ipoyo + e ipow^ p + + ^-ipow _ e ipo2/o^ p_ j ? (ni.46) 

which can be shown to approach to the continuum SF propagator [10] without any 0(a) 
term. 



D. Phase of Dirac determinant 

In general the determinant of the SF Dirac operator is not real for the overlap fermion 
since there is no 75 Hermiticity. Instead we have a following "Hermiticity" relation 

l5 uD SF uh 5 = 4 F , Ts^AWts = ^sf- (HI.47) 

However one cannot conclude reality from this relation since the SF Dirac operator connects 
different Hilbert sub-space as in (III. 21) and the determinant cannot be evaluated directly 
with -Dsf- We need to make a "Hermitian" Dirac operator which connects the same Hilbert 
sub-space in order to define the Dirac determinant. This is accomplished by 1^75 or w'75, 
which turns out to be 75 in the continuum. We define 

h sf = zWts = n +J D M t 75 n + : - w_, (m.48) 

H'sp = D SF u' l5 = U + Du' l5 U + : H„ -> H-. (111.49) 
The determinant is evaluated on the sub-space 7i- 

det H SF = det (u + Du^ 5 U + + U_) , (111.50) 

{H-} ^ ' 

det i^sF = det (U + Du'-f 5 U + + n_) , (111.51) 

{H-} 

where the right hand side is understood to be evaluated in the full Hilbert space by filling 
the opposite sub-space H + with unity. 

The phase of the determinant is given as follows 

( dfit } = e- 2l0W (det } H^j , (111.52) 
e- 2 ^ = det ( 75 m) 2 = detu, (111.53) 
e" 2 ^' = det (75^) 2 = det u ] = e 2 ^, (111.54) 
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which is not real in general. In the second equality of (III. 53) and (III. 54) we used a relation 



£ = ccdV, r = u; f Eu;, uj = e^ 75 . (111.55) 

The determinant of the unitary operator u is given by a product of eigenvalues \ n of the 
overlap Dirac operator 

detw= [] (l-oA n ), (111.56) 

ne{+} 

where product is taken over a sub-space in which the eigenvalue of E = +1 and the conjugate 
eigenvalue A* does not necessarily belongs to this sub-space. 

However we notice that this complexity of the Dirac determinant is not an essential 
problem since the phase is an 0(a) irrelevant effect and disappears in the continuum limit. 
Furthermore if we consider variation of the phase 

% % r / \ " 

5 e (x)<f> = 2 tr ^W MM_1 = ~4 atr [ T ' 5 e(x)D (l - aD ] j\ (III. 57) 

under a local variation of the link variable 

Se(x)Un(x) = ae^(x)U^(x) (111.58) 

we can show that S e ( x )4> is localized at the boundary. Since £ contains time reflection R and 
both of the operator 5D and (l — aD^ are local, the trace in (III. 57) has a contribution 
only at the boundary. Contribution from the bulk is suppressed exponentially by the locality 
property. 

For practical application to numerical simulation this phase problem should be settled. 
This is possible for even number of flavours by using two definitions of "Hermitian" Dirac 
operator (III. 48) (III. 49) and a fact that the phase is opposite for these definitions as shown 
in (III. 53) (III. 54). The phase can be absorbed into re-definition of the fermion fields. 
Explicit form of the two flavour Hermitian Dirac operator is 

H^=( DsF _ W75T 1 ' 2 : H-®H-^H-®H-, (111.59) 

where the unitary matrix is defined as 

U={ U u „) (111.60) 

to act on flavour space and r a is the Pauli matrix. 75 is also understood as two by two on 
flavour space. The Dirac operator is exactly Hermite for this definition: = H$ and 

the determinant is real. 

For single flavour case it is not still clear how to solve this practical problem. However 
the determinant is real at tree level and one may expect that the phase is settled as one 
approaches to the continuum limit. 
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E. Surface term 



When extracting the renormalization factors of fermions it is convenient to consider a 
operator involving the boundary source fields 



C(f): 

C'(x) 



Sp(x) ' 
5 

: 5p'{x) ' 



5p(f)' 
5 



(111.61) 
(111.62) 



where p, ■ ■ ■ ,p' are boundary values of the fermion fields given in (1.3) and (1.4). Coupling 
of the boundary value to the bulk dynamical fields was naturally introduced in the Wilson 
fermion [7]. However this is not the case for our construction since the boundary value 
vanishes with the orbifolding projection. 

In this paper we regard the boundary value as an external source field and introduce its 
coupling with the bulk fields according to the criteria: the coupling terms (surface terms) 
are local and reproduce the same form of the correlation function between the boundary 
fields in the continuum limit. Here we define the boundary values on the physical quark 
fields 



P+q(x)\x =o = p(x), P-q(x)\ X0=NT = p'(x), 

p'(x). 



g(x)P_ 



\xq=Nt 



xq=0 



L=o = p(x), q(x)P. 
One of candidates of the surface term is 

^surface = 3 ^(- p(x)P_q(x) \ Xq=q - q(x)P + p(x 

X 

- p'(x)P + q(x)\ xo=NT - q(x)P_p'(x)\ xo=NT ), 

where q and q are active dynamical fields on the boundary. 

According to Ref. [10] we introduce the generating functional 

Zf[p',p';P,P;V,V,U} = J D^Dijexp^-S F [u,ij,ij;p',p',p,p 

+a 4 ^] (il)(x)r](x) +f}(x)ijj(x) S j >, 



(111.63) 
(111.64) 



(111.65) 



(111.66) 



where r\ (x) and rj(x) are source fields for the fermion fields and the total action Sf is given as 
a sum of the bulk action (III. 19) and the surface term (III. 65). We notice that the fermion 
fields ip and ip obey the orbifolding condition (III. 9). We decompose the fermion fields into 
classical and quantum components 



*l>(x) = ipci(x) + x(x), if>(x) = ip cl (x) + x(x) 



(111.67) 



and insert it into the generating functional. The correlation functions between the boundary 
fields are derived with the same procedure as Ref. [10] by making use of this decomposition. 
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ip(x)ijj(y)) = G SF (x,y), 
71 



(q(x)q(y)) 

q(x)((y 
q{x)t{y) 
(((my)) 

(amy)) 
amy) 

at)C\y, 

amy)) = p 
am\y)) = p- 



1 - -D ) G SF 



1 - -D ) G SF 



-- P_ 



1 - 2^ ) Gsf 



+ 



1 - 2^ ) G sf 
1 - U -D ) G SF 



1 - 2^ ) Gsf 



1 - -L> ) G SF 



1 - -L> ) G SF 



x,y), 
(x,y)P+ 

(x,y)P. 

x,y) 
(x,y) 



J/o=0 
y =N T 
xo=0 



x =N T 

{x,y)P+ 



%v)P- 
[x,y)P+ 

(x,y)P- 



zo=0,i/o=0 
x =0,y =N T 
x =N T ,yo=0 



x =N T ,y =N T 



(111.68) 
(111.69) 

(111.70) 

(111.71) 

(111.72) 

(111.73) 

(111.74) 

(111.75) 

(111.76) 

(111.77) 



Here we adopted the physical quark field of (III. 15). The propagator Gsf is defined in 
(III. 37). We notice that the above propagators between the boundary fields and physical 
quark fields approach to the continuum SF boundary propagator without any 0(a) term. 



F. Mass term 

The mass term may be introduced with the same procedure as the continuum theory. 
We consider a mass matrix M which is consistent with the orbifolding symmetry 



TM + MT = 0. 



(111.78) 



Since the orbifolding transformation is the same as the continuum one on the physical quark 
fields, a naive candidate is to couple the continuum mass matrix mrj{xo) to the physical 
scalar density consisting of q(x) and q(x). Corresponding to various definition of the quark 
fields (III.15)-(III.17) we have several definitions of the mass term 



(111.79) 



where 77 is an anti-symmetric step function (11.14) on lattice. 

However we encounter a problem with this naive definition of mass term, since the 
massive Dirac operator does not satisfy the "Hermiticity" relation (III. 47). The phase of 
the Dirac determinant becomes mass dependent although it is still irrelevant 0(a) term. 

In order to avoid this unpleasant situation we propose even flavors formulation. For two 
flavors case we define the two by tow Dirac operator as 
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where 

Da^(m), = TT , f D 4- mrj (\ — ^ 



m(m)= MV n n , x , (IH.80) 
V D S F{m) 2 J 



D SF (m)i = n+ ( D + mr^ ( 1 - ) ) n (111.81) 



D SF (m) 2 = n + + m (l - |l>) ufyi 4 ) n_. (111.82) 
A "Hermiticity" relation can be found for this two flavors Dirac operator as 

4fW = r 1 ' 2 75 C/4 2 FV)^ t 75T 1 ' 2 , (111.83) 

where U is defined in (III. 60). 

The Hermitian Dirac operator can be defined to connect the same Hilbert sub-space as 

H§{m) = ^ 2 F ) (m)t /t 75T 1 ' 2 : H^@H^ ^H^®H-, (111.84) 

which is re-written in a trivially Hermitian form by a unitary matrix V 

The determinant of this Dirac operator is evaluated in a single Hilbert sub-space 

det HB(m) 

and becomes real. 

For the other mass term using the physical quark fields of (III. 16) and (III. 17) we also 
have two flavors Dirac operator with 

= ( H + (/J + m^t)nJ' ( IIL86 ) 

which satisfy the same "Hermiticity" relation (III. 83). 

Here we notice that U(2) flavour symmetry is broken to U(l) v x Z7(l) 3 by mass term 
like chirally twisted Wilson fermion. However the symmetry is recovered in massless limit 
and m = simulation is possible for relatively small box size we expect that this flavour 
symmetry breaking is not a serious problem. 

G. 75 mass term 

We have another candidate of the mass term which is consistent with the orbifolding 
symmetry (III. 78). That is the 75 mass term. In order to be able to define the Hermitian 
Dirac operator we need two flavors also in this case. The orbifolded action becomes 

S = la 4 ^D SF (m 5 )^ (111.87) 
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with the SF Dirac operator in the two flavors space 

D SF (m 5 ) = n+ (d + 27 5 r 3 m 5 (l - ^D^j fi_. (111.88) 

The Hermite conjugate is expressed as 

D SF (m 5 y = r 1 ' 2 7 5 w J DsF(m 5 )w t 7 5 T 1 ' 2 = t 1 ' 2 ^ D^im^u'^r 1 ' 2 (111.89) 

and we can easily show that the phase of the determinant is mass independent in this case. 
Furthermore we can absorb the phase into the fields and define the Hermitian Dirac operator 
in the same manner 

H SF (m 5 ) = J D S F(m 5 )f/ t 75 r 1 ' 2 , (111.90) 

for which we have exact Hermiticity ifsF^s)^ = -^sf(^5) and real determinant. 

Up to now the Hermitian Dirac operator (III. 59) (III. 84) (III. 90) contains gauge depen- 
dent unitary matrix U\ which may be an obstacle for practical application. This problem 
may be solved by chirally twisting the fields and adopting different boundary condition in 
temporal direction. We consider following chiral rotation of the fields 

$ = e -^ Br y, i> = ^' e -^ 75r3 . (111.91) 

The action is re-written as 

S = la 4 J2^D S F(m)?P', (111.92) 
D SF (m) = n_D(m)n_, D(m) =D + m(l- ^D^j , (111.93) 



where the projection operator H± is defined by using the time reversal operator £ of (II. 3) 
and third component of the Pauli matrix 

1 + Tr 3 

n± = —j— . (in.94) 

We should notice that the new fields ip' and ip' obey a different orbifolding condition 

U+ip'(x) = 0, i>\x)n + = (111.95) 
and satisfy a different condition at the boundary 

P + ^'(x) _ =0, P-^'(x) = 0, (HI.96) 



XQ=0 



x =T 



xq=0 



ij(x)P + _ =0, ij(x)P. = (111.97) 



x =T 



with a projection 



P± = r = t l5l0 r 3 . (HI.98) 
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Since the boundary condition is completely different, this formalism may not offer the 
same SF renormalization scheme. However it may be possible to define a new renormalization 
scheme with this theory. This formulation gives an equivalent partition function at least for 
the massless case connected by the exact chiral symmetry. This system has a gap caused 
by the discretization of p in (III. 29) and it may be feasible to make use of the same good 
property of the SF formalism in defining the renormalization scale. 

It is clear that this Dirac operator gives real determinant with Hermiticity relation 

7J SF ( m )t = 75 r 1 ' 2 D SF (m) 75 r 1 ' 2 . (111.99) 

The flavour symmetry is broken to U(l)v x ^(1)3 as other massive theories and SU(2)f 
symmetry in massless limit is guaranteed in a form of "Ginsparg- Wilson relation" 

75T 1 ' 2 ^ S f(0) +^ S f(0)7 5 t 1 ' 2 = aD SF (0) 75 r 1 ' 2 ^sF(0). (III.100) 

We notice that the projector Sr 3 commute with the Dirac operator 

"D(m),ET 3 ] =0 (III.101) 

and the eigenvalue of this new system becomes the same as that of the ordinary massive 
overlap Dirac operator D{m) with half numbers of degeneracy. The propagator is defined 

as 

**■,(*,„)= 2 (D*(m)-) = 2 (n_-i- T n_) , (III.102) 



x,y \ Dim) 

\ v / / x ,y 



which takes a simple form at tree level 



g laN T n= _jv T+ i + M 2 



Xe i P0 xo I ^i Po y + e i Po yo^ p + (e-ipoyo _ jPoyo^j p + j ? (HI. 103) 



where 



A„ = (l - ^rnj Ap, M = m+(l- ^rnj B (III. 104) 

and A M , B are defined in (III. 35). The temporal momentum p has a discretized value of 
(111.29). 



IV. CONCLUSION 

In this paper we propose a new procedure to introduce the SF Dirichlet boundary con- 
dition for general fermion fields. In the former formulation the boundary condition is in- 
troduced by cutting the Dirac operator at the boundary [7]. For the Wilson fermion the 
boundary condition is automatically decided depending on the signature of the Wilson pa- 
rameter r. However this formulation produces an exponentially small eigenvalue in the kernel 
D w of the overlap Dirac operator because the relative signature of the Wilson parameter 
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r and the mass parameter M is opposite. This may cause breaking of the locality of the 
overlap Dirac operator or breaking of the axial Ward-Takahashi identity in the domain-wall 
fermion. 

Instead of cutting the Dirac operator we focus on a fact that the chiral symmetry is 
broken explicitly in the SF formalism by the boundary condition and adopt it as a criterion 
of the procedure. We also notice that an orbifolded field theory is equivalent to a field theory 
with some specific boundary condition. We search for the orbifolding symmetry which is not 
consistent with the chiral symmetry and reproduces the SF Dirichlet boundary condition on 
the fixed points. We found that the orbifolding S 1 /Z 2 in temporal direction including the 
time reflection, the chiral rotation and the anti-periodicity serves this purpose well. 

Application of this procedure to the overlap Dirac operator is straightforward since this 
system has both the time reflection and the chiral symmetry. The eigenvalue of the Dirac 
operator and the propagator are derived at tree level, which are shown to agree with the 
continuum results in a — > limit. 

We found a technical problem that the Dirac determinant is complex. However this is 
not essential since the phase of the determinant is an irrelevant 0(a) term and its variation 
is localized at the boundary. Furthermore we can absorb this phase into re-definition of the 
fermion fields for even flavours case and explicit form of Hermitian Dirac operator was given 
for two flavours. For single flavour case it is not still clear how to deal with this problem. 
However the determinant becomes real at tree level and one may expect that the phase is 
settled as one approaches to the continuum limit. Numerical simulation is needed to observe 
scaling behavior of the complex phase. 

The mass term may still have a problem. If we define the mass term in an ordinary way to 
couple to the physical quark field with kink function t](xq) the phase of the Dirac determinant 
becomes mass dependent. In order to avoid this problem we need even number of flavours, 
for which the phase can be absorbed into the fermion field. The chirally twisted mass term 
serves for the same purpose. This mass term is consistent with the orbifolding symmetry and 
gives mass independent phase of the determinant, which can also be absorbed by re-definition 
of fields. The most satisfactory solution may be to introduce the chirally twisted boundary 
condition, for which the Dirac determinant is naturally real without any re-definition of 
fields. We should notice that the U{2) flavor symmetry is broken to U(l)v x U(l)f for these 
two flavors formulation. However the symmetry is recovered in massless limit and m = 
simulation is possible for relatively small box size the flavour symmetry breaking may not 
be a serious problem. 

The SF formalism with the domain-wall fermion can be formulated in the same way 
[35] since the symmetry is exactly the same as the overlap Dirac operator [28]. Application 
of this procedure to the Wilson fermion is not trivial since the Wilson term does not anti- 
commute with T operator but commute with it. In order to be consistent with the orbifolding 
symmetry we may need to introduce the step function rj of (11.14) in the Wilson term. 
However the Wilson term vanishes on the boundary x = 0, T and the doublers appear. 
We may need to vanish all the fields at the boundary, with which we reproduce the same 
SF Dirac operator as the original one [7]. Introduction of 75 into the Wilson term may be 
another way to apply this procedure. 
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